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SUMMARY 

An axial -discharge mixed-flow compressor, which is especially 
adapted for jet engines "because of the large mass flow per unit 
frontal area, is described. The basic concept of the design pro- 
cedure is the incorporation of conditions for efficient flow into 
the design equations so that all reasonable cases can be computed 
without investigating the effect of a large number of design param- 
eters. General equations of relative fluid motion are developed to 
show clearly the assumptions involved, and the empirical character 
of the simplifications employed to render the system of equations 
solvable and to complete the design. 

The best impeller was selected on the basis of the maximum 
air-flow capacity, which was 19.6 pounds per second for a 14-inch- 
diameter impeller with a tip speed of 1480 feet per second and a 
pressure ratio of 3.5. 


INTRODUCTION 

The centrifugal compressor has the advantage of high pressure 
ratio in a single stage, simple construction, mechanical strength, 
compactness, and reliability, but is handicapped by relatively low 
efficiency and low mass flow per unit frontal area. For high-speed 
jet engines, the frontal area is important from the standpoint of 
drag and mass flow is important from the standpoint of power. The 
axial -flow compressor, on the other hand, has the advantage of high 
efficiency, high air-flow capacity per unit frontal area, and ease 
of staging, but is complicated and expensive to manufacture, frag- 
ile, and relatively heavy. A compressor combining the best features 
of these two types would be rendered compact, simple, light, and 
strong by accomplishing the compression in rotors of high blade 
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solidity and work output, and would be reduced to a small frontal 
area by having the air directed axially at the entrance and exit 
of the impeller and in the diffuser. The apparent lower efficiency 
of the centrifugal compressor could probably be improved by the 
development of a more rational method of design, thus Justifying 
the development of a compressor with high blade solidity and com- 
pression in a single stage. 

A compressor of this type was therefore developed with a 
maximum internal diameter of 14 inches including the diffuser, an 
equivalent tip speed of 1480 feet per second based on stagnation 
inlet conditions, and a compression ratio of 3.5. The air flow 
was to be maximized in the design process. 


AERODYNAMIC BASIS OF DESIGN 

The basic problem of the design is to ensign the correct, 
velocity distribution at the impeller entrance and exit and to 
select dimensions that will produce maximum weight flow for a 
prescribed outer diameter, rotative speed, and pressure ratio. 

A study of the effect of all design variables on air -flow capacity, 
however, is a stupendous task because of the large number of cases 
that must be computed. Relations between the design parameters 
are therefore established on the basis of certain assumptions as 
to how the air could be efficiently handled and the independent 
variables are reduced to two, thus permitting an evaluation of all 
reasonable designs with minimum effort. 

Equations of flow at the entrance and exit are first developed 
and these flows are related by means of conditions for efficiently 
handling the air. This part of the design is completed with a 
maximized flow before the Impeller shape is determined. An 
impeller shape Is then computed to accomplish the desired change 
in air flow and state from the entrance to the exit. General 
equations of motion in the impeller are developed In order to 
Indicate clearly the basic assumptions made in obtaining the shape 
of the impeller hub. Appropriate simplifications are used to per- 
mit an easy solution for the impeller design. St re am -filament 
methods are used in determining the velocity distribution and 
blade shape in the inducer section, and from this point downstream 
a faired curve Is used for the blade camber line. 
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Basic Assumptions 

The following assumptions are used in establishing flow rela- 
tions and velocity distribution; 

1. At the entrance and exit (stations 1 and 3 in fig. l), the 
entropy is constant on an axial plane, which is the equivalent of 
assuming no boundary layer on the shrouds and equal entropy increases 
for all streamlines in the impeller. 

2. The gas inside the impeller flows on surfaces of revolution. 

3. Tangential velocity at the casing in the discharge annulus 
is 0.95 of the tip speed. 

4. The velocity reduction at the blade tips in the inducer 
section is one -third the incoming relative velocity. 

5. Pressure rise allowed on the case is slightly less than the 
centrifugal pressure rise on the assumption that the boundary layer 
limits the pressure rise. 

6. At the exit, the maximum allowed absolute Mach number 
is 1.4. 

7. At the entrance, the maximum allowed relative Mach number 
is 1.0. 


For computing the shape of the impeller there are the follow- 
ing assumptions ; 

1. The blade tips are shaped according to two-dimensional 
stream-f i lament theory, which is applied to compute the velocity 
on the blade surface. Suitability of the velocity distribution is 
judged by two-dimensional boundary -layer theory. The rest of the 
blade, extending down to the root, is shaped to maintain radial 
blade elements. 

2. Inside the Impeller the entropy is constant on surfaces of 
revolution normal to the meridional projection of the streamlines. 

3 . The frictional forces on only the blades and impeller are 
taken into account; those on the case are neglected. 

4. An infinite number of blades is assumed for computing the 
hub shape. 



4 


MCA RM No. EOF 04 


Determination of Entrance and Exit Conditions 


At the entrance, it is assumed that radial equilibrium exists 
at the impeller face. For balance of the centrifugal and pressure 
forces 


where 


1 dp _ V u 2 

p dr r 


(X) 


p pressure 

r radius 

V absolute gas velocity 
p density 

Subscript : 

u rotational component 

From thermodynamics, with the assumption of constant entropy, 

dH = I dp (2) 

P 

and 

V 2 

H - - V (3) 


where 

E enthalpy 

Subscript : 

t stagnation state 

Substitution gives 

dH _ _ _d_ /_v£\ _ V 

dr “ dr dr \ 2 / 


r 


( 4 ) 
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At the entrance, the, stagnation enthalpy is constant because of the 
constant energy level, and (4) reduces to 



JL (rV 

2r 2 dr ( 


( 5 ) 


where the subscript a indicates the axial velocity component. The 
mass flow W between the cylinder of radius r and the case r c 
is 


P V a rdr (6) 


The impellers considered are to have subsonic internal relative 
flow and the relative velocity at the entrance blade tips is there- 
fore assumed to be sonic. The mathematical expression for this 
condition is 





2 

7+1 


( at ,!) 2 + 



■ (^U,c,l r c,l) /( r c,l) 

where 


( 7 ) 


a sonic velocity 

7 ratio of specific heats 

(jo angular velocity 

Subscripts : 
c case 

1 impeller entrance 

Station numbers and some dimensions and flow parameters are shown 
in figure 1. For convenience, the symbols used in the text and 
appendix A are listed in appendix B. 
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With these three equations for the entrance flow, insufficient 
relations are available to establish the flow there. At the exit, 
additional conditions are required to prevent the work output and 
exit Mach number from varying too much. In order to apply equa- 
tions (5), (6), and (7) it Is therefore necessary to relate the 
entrance to the exit flow and to apply the conditions imposed on 
the exit flow to determine the entrance flow. 


Flow on Impeller Blade Tips 


The gas near the blade tips is assumed to be compressed to 
such a degree in the inducer section that the boundary layer on 
the blade tips will separate from the blades with a further pres- 
sure rise in a reasonable distance. Downstream of the end of the 
Inducer, a pressure rise in the radial direction slightly less than 
that due to centrifugal force is therefore allowed at the case. 

In symbolic notation, 

dp < poo 2 r c dr c 

There is an entropy rise of the air and the pressure rise is 
therefore arbitrarily reduced 


1 dp 2 
- = oo r 

p dr c 



where 

S entropy 

T absolute temperature 

But 


1 dp T d S + dH 
p dr c " ' dr 0 dr 0 

and Euler’s equation combined vith the energy equation gives for 
any streamline 

dH-k = ood(rV u ) = d^H + " 2 “^) 
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Therefore, 

2 “ d ( r c v u,c) d fro 2 ] 

" r ° = dr 0 ' dr c V 2 / 

If the relative velocity v is introduced, the solution of the 
equations is found to he 


t c = constant = v Cj2 = v Cj3 


( 8 ) 


vhere the subscript 2 indicates the end of the inducer section 
before radial flow begins and 3 indicates conditions at the 
impeller exit. A compression in the inducer is estimated as 
f easible with 


c , 2 = 


2 

5 


c,l 


(^a,c,l) + (^u,c,l " “ r c,l) 


„ll/2 


( 9 ) 


Flow at Impeller Exit 

If the blade tips at the exit are directed axially at this 
point the tangential component of the velocity is 


^u,c,3 “ v u,c,3 + wr c,3 ~ * wr c,3 (10) 

where f is the slip factor. For the impellers under considera- 
tion, eighteen blades are used, and an estimate gives f = 0.95. 

The axial-flow component is then 


a,c,3 


= v 


a,c ,3 = 



v u,c,3) 



- (1-f) 2 to 2 (r c>3 ) 2 


( 11 ) 


In the discharge annulus the velocity components will vary in some 
manner from the values at the tip. If a linear variation of the 
quantity (V u ^3) 2 /r 3 is prescribed, 


(Vs) 2 


(Vh. 3) 2 


v.-j- 

r h,3 


l^u,c,5^ / r c,3] ~ 


c,3 


“ %3 


(r-7 - 


h,3 


( 12 ) 
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where subscript h indicates the inner -radius (hub) or blade -root 
value. The enthalpy at the exit is a fundamental compressor per- 
formance parameter, and is, of course, related to the exit veloc- 
ities. Because of limits to be imposed on the exit gas velocity at 
the blade roots, the enthalpy at the exit will be related to the 
Mach numbers of the absolute and relative flow at the blade roots 
in the discharge annulus. If equation (4) is now applied at the 
exit and (12) is substituted, the resultant equation can be inte- 
grated to give 


h 3 " E h,3 = 


or 


^u,h,3 / v2 u,o, 3 V 2 u,h,5 \ ( r 3 ~ 711,5) 

r h,3 \ r c,3 ' *h,3 ) 2 ( r c,3 r %3 


-,3] 


(13) 


E c, 3 " E h,3 


r c ,3 " rh,3 


( v u,h,3) (^u,c,3) 

2 


L %3 ' %3~ J 


(14) 


The Mach number at the exit must not be too high for efficient 
diffusion. Its value is given by 


„ 2 V 3 2 V // J^_\ 

”3 - = (7-l)H t)3 /[ 2H tj3 y 


(15) 


The relative Mach number M ? is also of significance because its 
magnitude will indicate choking flow in the impeller. It is 
related to other flow variables by 


(M's) 


2 


(v 5 ) 2 _ (Vs) 2 + ( y u,3 - M'S) 2 

(V-i) H 3 = (r-D % 


(m- 3 ) 


2 



“ 2 (r 3 ) 2 - 2 w 3 V Uj3 

(7-1) h 3 


(l^T 
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The value of H ?j at the root (H^ 3 ) from (16) is substituted in 

(14) to give the value at the case related to the Mach numbers at 
the root 


H, 


c,3 


2w h,3 Vh,5 - ^(r h|5 ) 2 


(7- 1 ) (M h>3 ) 2 - (M' hj3 )2J 
(r, 


- r h,5) | /(Y u,h,5 ) 


V 


h,3 


(V ) 2 
r c,3 


(17) 


There are now available sufficient relations to determine the 
complete flow at the entrance and exit with only two independent 
parameters, provided there are reasonable limitations on the flow 
conditions. A number of designs are therefore computed with 
various assigned values of these parameters and the design is 
selected on the basis of the best air-flow capacity. The assump- 
tions and procedures for computing the entrance and exit velocities 
are now summarized. 


SCHEDULE OF DESIGN COMPUTATIONS 

For all impellers to be designed, the work output per pound 
is assigned and known. Consequently, from the state of the enter- 
ing gas (p t ,, Hj. j) , c 3 is known. The work output varies 

from streamline to streamline, and therefore the mean work output 
per pound is known only approximately. The rotative velocity go 
and the maximum frontal dimension r c ^3 are also assigned. 

Computations at Case 

1. From the assigned data, v u c 3 and V u are computed 

by means of ( 10 ), which assumed axially directed blades at the 
exit near the case. The Euler equation 

Ht,3 ' H t,l = “( r 3 V u,3 ' r l T u,l) US) 

# 

can then be used to find the entering moment of momentum at the 
( r c, l v u,c, l) * 


case 
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2. The particular design being investigated now enters into 

consideration by assigning values to r c,l r h, 3 * 

assigned value of r c p permits the computation of V u 0 ^i from 
(r c ]V U c p) . The value of r h ^ 3 is only tentative because in 

some cases the entrance annulus will not accommodate the entire 
exit flow. Therefore, when the minimum permissible r h p is 

computed, the corresponding streamline position at the exit must 
be found and the exit annulus cut off at that position. 

3. Equation (7) is used to find V a c p. A relative Mach 
number of 1.0 at the entrance blade tips is assumed. 

4. For the exit at the case, equations ( 8 ) and (9) give v c ^p, 

v o, and v„ ,. With v - and v.. , equation (11) deter - 

C, C. J C,0 0,0 U.,Ay,vJ 

mines Y 0 „ , which with V„ r 3 is used to find V c 3 . Then 
a, c,o u.,o,j , 

the enthalpy H c ^3 is computed from 

H = H t - tL (19) 


By assigning an efficiency for the impeller (0.85), it is possible 
to find p t c 3 and P c ,3 Prom 


H 


\7-l 


Pc, 3 


c .3 V p t,c,5 




RT t ,c,3 


p t , c , 3 1 H ' 


C-L 3 


1 

7-1 


7-1 H t , c , 3 \ H t,c,3 / 


( 20 ) 


where R is the gas constant. 


Computations at Exit 

5. The first step in computing the flow conditions and gas 
state at the exit is to find the conditions at the blade roots . 
Because the Mach number increases with decreasing r 3 , and because 
the efficiency of normal shock compression drops rapidly for Mach 
numbers increasing above 1.4, a limit of 3 = 1.4 is taken at_ 

the root. There are two possible procedures for the next step. 

(a) The maximum flow at the blade roots is assumed by 
setting M'k 3 = 1.0. With these values for the relative and 

absolute Mach numbers, equation (17) then determines v u ^h,3> 
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and % 3 is then computed from equation (16). Equation (15) 
is then employed to find V^3 and 3; is com- 

puted by a right triangle formula. In completing a design begun 
on such a basis, it was found that the flow limitation occurred 
at the entrance rather than the exit and the condition 
M'k 3 = 1.0 was therefore discarded. 

(b) An alternative procedure assigns several values of 
Hfc Yi 3 for ©ach r h 3. When the entrance flow is computed for 

each case, it is found that this flow could be accommodated at 
some radius for the root at the exit. Of the values of the 
exit root radius thus found, that which is equal to the assumed 
r^ 3 then gives the desired Ht^h,3« For any assumed value 

of k 3, the procedure would consist of first finding 
and 3 from equation (15). Equation (14) is solved for 
V u h 3 and V & h 3 ^hen found, thus determining all needed 
flow conditions at the blade roots. 

6. In computing the variation in fluid state and flow condi- 
tions at the exit, 3 is assumed to vary linearly with radius 

between the value at the case and that at the blade roots. Equa- 
tion (13) gives H3 and (12) determines V u ^3. Then 


(v 3 ) 2 = 2(11^3 - H 3 ) 

(V 3 ) 2 = (V 3 )2 - (V u,3 )2 

and by assuming isentropic conditions at the exit 

“5 ‘ (£)' 1 < 21 > 

For correlation of the entrance and exit streamlines, the mass flow 
V r 3 between the case and a cylinder of radius r^ is required. 



W = 2ir 


p 3 v a,3 r 3 ^3 


( 22 ) 
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Computations at Entrance 

7. Computations of flow at the entrance use equation (18), 
which gives 

r l v u,l = S ^,1 “ + r 3 V u,3 ( 23 ) 

The entering moment of momentum r^ Y u ^ is therefore known as a 
function of r 3 , the radius at which that element of gas leaves 
the impeller. However, the value of the mass flow between 

radius r-^ and r c ^ is the same as that between r 3 and r Q ^ 
for the same streamline; that is, 

1 

p l V a,l r l te l 

1 
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d ( r l Y u,l) 2 ' _ 


dr 


3 d 


dW 


- 2(r lV i) dw d I r 3 V u,3 




CO 


But 


dr. 


dW 2*r 3 P 3 V a;3 

2 

Also, 3 and (V u 3 ) A 3 are linear functions of r 3 ; 

sequently 

^,3 _ ^t ,5c ~ ^t ,3h 


con- 


dr. 


•0,3"% 3 


and 


d^ <^«,3> = ^ 


. ( r s) 2 R Y u,c ,5^/ r c , 3 ] - [(Vg.h.s) 2 /^,! 

2 Y u,3 r c,3 " ^,3 


Therefore 


d (n ? u.i) 2 _ nVi 

M " dW nr 3 P 3 V a,3 


x< 


,K,3 

r 2 
3 

^u,c,3> 2 

(v ) 2 

v u,h,3 ; 

- ^ i c |3 th,^! 

^ 2 

2 ^u,3( r c,3 “ 3^ 

r -z 

L c,3 

%3 

( r c,3 - r h,3)j' 


(26) 

can be evaluated from exit computations as a function of W before 
solving for the entrance flow. Now equation (25) can be written 


d ( V a,l> : 


aw 


( r i ) 


27 ) 


and equation (24) can be expressed 
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a(r x ) 2 
dW = 

The system of equations (27) and 
by the Kutta-Runge method. 


* p l v a,l 


(28) 


(28) can he numerically integrated 


Another scheme of good accuracy for this special system was 
devised. Equation (2?) can he expressed 


d ( y a.l) 2 _ 1 

d ( r l v u,l) 2 < r l) 2 


where for convenience 


d 2 y 1 

. 2 , ,4 

dx (r x ) 


dy 1 

dx - (ri )2 

y = ( V a,l )2j X = (r l V u,l )2 - Th8n 

d ( r i ) 2 i aw d ( r i ) 2 


(r x ) 4 d (r x V Ujl ) 2 dW 


£jL m i 

dx 2 np l v a ld( r l) 4 

A Taylor' 8 aeries expansion for y is 

y (x+A ) = y(x) + Ay’(x) + y M (x) + ~ r y ,M (x) + 

u O • 

or 

y(x-A) = y(x) - Ay' (x) + £ y"(x) - y’" (x) + + 


(29) 
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Then 


and 


y ( x+A ) = 2y(x) - y(x-A) + A 2 y"(x) + + 


y(x+A) = y(x-A) + 2Ay ' (x) + ~ A 3 y"* (x) + 


(30) 


(31) 


2 

In starting the solution, a small increment in (r-jV u 1^ use< ^- 

end the term ^ A° y"' M i n equation (31) is neglected to obtain 


(Va,l 2 ) b * (V a ,l 2 ) a = - 


(v u ,x 2 n 2 )a - (v u ,hnhtj 


2' 


Jn 2 ) a + ( r i ) b 


(32) 


as the difference in axial velocities squared at two inlet radial 
stations a and b for the estimated value of (r-^)b. The mean 

value of p-jV p -] for the step is then estimated, utilizing the 


guess for 


a? 1 

( r i^)b 


to compute the value of at the average 


position by means of the stagnation density and the velocity. With 
this mean value of p^V a ^ for the step from a to b, a better 

2- 


estimate of (r-,^)h is obtained from 


(n 2 ) b 


_ ( r 2\ 1 

1 a ' "Pl V a,l 


(Wv, - Wj 


(33) 


where V-^ and W a are known in advance for the corresponding step 
In (r£V u p)-b “ ( r i v u l^a* T ^e nev value ( r l 2 ^h can ^en be 

used again in equation (32) and the process repeated to find the 
value of (r-^ 2 )^ that checks. After one small step has been 

taken, equation (30) Is used because It neglects only the fourth 
and higher derivatives. The size of the Interval can be doubled 
after the second step because of the higher accuracy of the formula. 
In the case of equation (30), the steps are uniform In (r-jV u 1 ) 2 . 

Then equations (30 ) and (29) give 


( V a,l 2 )d - z < v a ,l 2 )b - (V a ,l 2) a + 


r 


x 


( r i 2 v u,i 2 )d - (n 2 v u,i 2 ) b 


n ^ p l v a,l pr l ^b 

2 (34) 
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and parabolic (second degree) integration of (28) gives 


< r i 2 )d - < r i 2 >b = - 


(n 2v u,i 2 )d - h 


12 n 


X 


^ p l^a,l^a ( pp l7a,lA + ^ p l V a,l^d 


(35) 

In equation (35) , (V & ^) d and are known, but (pp)^ is not. 

If an estimate is made, then (r^ 2 )^ an& hence 


<V> = ^ (V u,l r l) 

and finally p, can be computed. Equation (35) is repeatedly 
used until the estimated and computed values of (r-^J^ agree. 

The method fails if (r^V u ^) as a function of W is constant; 

in that case however, a much simpler method of solution Is avail- 
able because equation (25) indicates (V & ^)^ is constant, and 

(V u ^) 2 and hence p are known functions of r^. The solution 

of equation (28) by Simpson’s rule then Identifies the variation 
of r-^ with variations In W. 

8. At this stage of the computations the entire velocity dis- 
tribution at the entrance and the exit is known, and the mass flow 
is also known as a function of radial position at the entrance and 
the exit. Examination of the velocity diagrams indicates in some 
cases a radius at which it is desirable to set the blade root. 
Computation of the blade stress at the roots also Imposes limita- 
tions. If the blades are tapered with increasing thickness toward 
the root, blocking of the flow area may Indicate a desirable radius 
for the root. By such considerations, a blade-root radius is 
decided upon and a mass flow established for the assumed values 

of r and r, ,. -- 

c,l h,3 


DESIGN OF IMPEELER 

After the best gas velocity distribution has been computed at 
the entrance and exit, the next problem is to design an impeller 
between these planes to accomplish the desired turning of air. 
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Inducer Section 

The inducer section is regarded as that portion of the impeller 
where the air receives an initial compression with no radial flow 
at the case. At the end of the inducer section, the flow area is 
made as large as possible by designing for no radial flow at the 
root. Conditions at the inducer exit are estimated by the dif- 
ferential equations of motion on the assumption of a certain pre- 
scribed relative velocity at the case. (See equation (18).) 

Inside the impeller 


% = H + = E tfl + (i>(r 2 V u ^ 2 - r^V^x) 


The quantity 

E + iv 2 - CorV = H 4 v 2 - ~ w 2 r 2 = x - oor.V (36) 

2 u 2 2 u J- u j x 


is therefore constant for any one streamline, and the relative 
stagnation enthalpy E + (v 2 /2) changes on any streamline only 
as a result of change in the potential energy level (u 2 /2). 
Equation (4) is then 


an (v u ) 2 

dr r 



(n v u,i) 




(rV u ) 


1 av£ 

2 dr 


(37) 


In terms of relative velocity 




2 (TV. 


u 





dW 

dr 


(38) 


Continuity gives 


dW 

dr 


-2jtr p v 

a 


In accounting for the area blocked by the blades, the quantity 
2irr - (Bd/cos must be used for 2nr where 


B number of blades 


d 


blade thickness 
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\|/ angle between direction of axis of rotation and intersection 
of blade surface with circular cylinder about axis of 
rotation 

The corrected mass -flow equation is 

e - (zicr - -■- " d -7- X ) pv cos y\i (39) 

dr \ cos j 

For structural strength it is prescribed that radial lines through 
the blade tips be blade elements and consequently 

tan \J/ - (r/r c ) tan \J/ C (40) 


Isentropic conditions are assumed at any axial depth; from this 
assumption the density variation is estimated by means of 


/ H ' 

p = Pc fc 


7-1 


in conjunction with (36) rewritten 


1 7 (jo^r^ 

H = Sfc,! ' wr l v u,l " 2 T + ~2~ 


(41) 


(42) 


The equation of motion is then 


a(y g /g) = 

dvr 


(2itr - Bd sec \j/) pv cos ^ 


sin' 2 \|/ + 2v oo sin \|/ 


r d < r i T «,i) 

- < 43) 


This equation may be solved simultaneously with equation (39) by a 
step-by-step process such as the Kutta-Runge method. The quantity 

(j o — (r-jV.. -i) is a known function of W and is computed before 
dw 1 U .? J - 

the inducer section is designed. The boundary value for v 0 = v r 0 
is also known, but \|/ is not. A solution is therefore obtained 
for several values of '\s n and the one giving the desired value for 
Is selected. 
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The blades must be turned from the entrance direction rapidly , 
especially in the beginning because the blade thickness may cause 
a velocity rise before the diffusion process is begun. The veloc- 
ity distribution at the blade tips is estimated on an assumed 
camber line by the stream -filament method as outlined in refer- 
ence 1. A check for boundary -layer separation is then made from 
the calculated velocity distribution and the suitability of the 
assumed design so evaluated. Modifications in shape are indicated 
by the desired change in velocity. The rest of the inducer blade 
surface is shaped according to equation (40) . 


Mixed-Flow Section 


From the end of the inducer section to the impeller exit, a 
smooth surface of revolution is assumed for the shape of the case, 
and on this surface the curve representing the blade shape at the 
tips is assumed. The blade surface is then obtained from equa- 
tion (40). The quantities r c and \J/ c are then known functions 

of the axial depth, and therefore \[/ is known at every point. The 
boundary condition of constant v c is attained by properly shaping 

the root. The equations utilized in this step are derived in 
appendix A, and are listed here for continuity. The mass flow 
between two surfaces that are meridional projections of the stream- 
lines is 

- dW = pv ( 2jtr - -22; — ) cos cl dn (44) 

\ cos a I v 


and the variation in relative velocity is determined by 


where 


V 


di(v c /2 ) 


v 


dW P (2„r - -JSJ) 

\ cos a j 


cos a | 

(2 co . sin ql\ 

1 sin a cos cp 



\v + r ) 

cos a 

dW 


( 45 ) 


a angle between velocity and meridional component 

dn element of length in meridional plane normal to meridional 
velocity component 
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radius of curvature of meridional projection of streamline 
qp angle "between meridional component and axis of rotation 
If the blade consists of radial elements then 


tan y\i = tan ^ (46) 

r c cos cp 

c 

Fluid density is related to that at the case and the enthalpy of 
the fluid by 



(47) 


The density at the case p is computed on the assumption of 

c 

constant efficiency during the compression process. Equation (46) 
is evaluated on a plane normal to the axis, and equation (47) is 
evaluated along the curve n, which consists of the line 
elements dn. An infinite number of blades and constant entropy 
along the n -curves are assumed. This system of equations is 
solved by taking small increments in the mass flow and evaluating 
the position of the streamline next to and inward from the one 
already known. The velocity is also evaluated at this streamline. 
The new streamline position is computed throughout the whole 
impeller before moving on to the next streamline. By this method 
it is possible to evaluate Rj^ and solve the equation. When the 

entire mass flow is thus accounted for, the shape of the hub is 
outlined by the last streamline. 


Impeller Design Dimensions and Performance 

The foregoing system of computations resulted In an impeller 
with the following design dimensions and performance character- 
istics : 
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Dimensions : 


Radius of case at entrance, inches 5.62 

Radius of case at exit, inches 7.00 

Radial blade clearance, inch 040 

Radius of root at entrance, inches 2.75 

Radius of root at exit, inches 5.90 


Flow Characteristics at Exit ; 


Absolute Mach number at case 1.28 

Absolute Mach number at root 1.40 

Relative Mach number at case 60 

Relative Mach number at root 75 

Rotation velocity component, feet per second 1400 

Axial velocity component at case, feet per second 710 

Axial velocity component at root, feet per second 840 


Impeller velocity at 7 -inch radius, feet per second 1480 

Perf ormance : 


Average pressure ratio 3.5 

Mass flow, pounds per second 19.7 

Efficiency near case 85 


The conditions at the exit are shown in figure 2 and those at the 
entrance in figure 3. Because of the constant entropy assumed at 
the exit the efficiency of compression is lower for streamlines 
near the root because of the lower work input. 


EVALUATION OF DESIGN PROCEDURE 

The fundamental objective of the design method was to find 
sufficient reasonable restricting conditions on the types of flow 
so that the number of variable design parameters would be as small 
as possible. The highest air-flow capacity design for variation 
of these parameters could then be found with a reasonable amount 
of labor. Seme of these restrictions are quite arbitrary, such as 
the energy conversion possible in the inducer and the limitation 
of the absolute Mach number to 1.4 which was chosen because of the 
high efficiency of normal shock compression at this value. How- 
ever, efficient diffusers can be designed at higher Mach numbers by 
utilizing oblique shocks. 

More exact calculations of flow detail are also desirable, 
especially at the impeller entrance where radial flows exist. 
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More detail of the velocity distribution on the blade surfaces inside 
the impeller is required to account for three-dimensional flows. To 
utilize such a knowledge to the utmost, however, solution of the 
equations for the boundary -layer flow on the rotating blades is 
required in order to decide what is a good velocity distribution. 

No such solution exists. 

One can therefore say that the design procedure is partly 
rational and partly empirical. Because the analytical solution of 
these flow problems is not in immediate prospect, experimental 
techniques must provide some of the information required for the 
step-by-step improvement of this promising type of compressor. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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APPENDIX A 


EQUATIONS OF INTERNAL RELATIVE FLOW IN IMPELLER 

The gas inside the impeller is assumed to flow on surfaces of 
revolution, which permits an extensive simplification of the equa- __ 
tions of motion. For_a compressible fluid with absolute velocity V 
and steady velocity v relative to the impeller, which rotates with 
angular velocity a), the equation of motion is 

2 _ 

vx(vxv) + ^ = (I ' 2 r + 2v x (J) + F/p (l) 


The bar over the symbol indicates a vector quantity; F is the 
force per unit volume exerted by viscous forces and bodies dis- 
tributed in the fluid. The tern 


2 ~ 

a) 41 r = 


u x 0) = or 




is the centrifugal force where 


u = co x r 
V = u + v 
2v x (0 


linear velocity of rotation 
absolute velocity 
Coriolis force 


From the laws of thermodynamics, the equation will be modified by 
expressing the pressure and density in terms of other gas -state 
functions 


i vp = vH - TvS = 

P 

In the section on inducer design, 
h = H t - V u 0)r * 2 t 


Ht - 


- TVS 


it was found that the quantity 


- = v 2 u 2 

- u*Y - H + £- 


( 2 ) 


was constant on streamlines in the inducer. This quantity is 
Introduced into the equations of motion in order to determine the 
conditions for which It is constant for streamlines In the impeller. 
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(l/p) vp = vh - TVS + (vu 2 /2) - vv 2 /2 
Equation (l) then hecoines 

vh = TVS + v X (2oo + v x v) + F/p (3) 

This equation shows that the condition for h constant on a 
streamline (v*vh = 0) is 

T al - - F v/p (4) 

where b/ba is the derivative along a streamline and F v is the 

component of F parallel to v. The entropy rise is a result of 
friction, and hence F v is a frictional component of force. The 

rest of the force is normal to the velocity vector v and is 
designated as a lift force Fj. Because h is constant along any 

streamline, then changes in h may he determined from values on 
different streamlines at the impeller entrance. At the entrance h 
is constant for different streamlines at the same radius. Because 
the flow is assumed to he on surfaces of revolution, h is there- 
fore constant on surfaces of revolution; that is, 

u*vh = 0 

The only remaining component of vh is therefore normal to v 
and u and is parallel to the unit vector 

A . V XU . (5) 

/(v X u) £ 

which is in the meridional plane and normal to the projection of v 
on that plane. If the assumption is made that vs is parallel to 

Tvs = - (%F T )/p 

and 

F z /p + v x (2co + vxv) = n|“ = vh 
where ^ is a unit vector parallel to v. 
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It Is now assumed that the lift force Fj is on the surface 

of revolution formed by the stream sheet as well as nomal to the 
velocity. Thus, F i is normal to ft, and therefore, 

^ = ft. y x ( 2(0 + v x \) ( 6 ) 

on 


Equation (6) provides the basis for computation of the streamlines 
in the impeller. 

For greater convenience the velocity is broken down into two 
components v m the meridional component and v u the rotational 

c omponent . Then 


ft-v x go = v*(o x ft = v • u cos cp/r = v u u cos cp/r = v u cocos cp 


Also, 


and 


v x v x v ss v v^/2 - v*vv 


VW = (v m + v u ) • V (v m ft + v u ft) 




2 c)ft / \ iu . 

= v w r — -f m v_ -r — + + etc. 


a , r m 


v^“ _ dv. 


A 'm — A ^ Y u a v u d y m 

sn ““ + v n¥ +ilT ii l ¥ + r ¥ 


R 


m 


v_ _ 
m 


ft 2 u 


6 dv. 


u 


+ ~ y u sin cp“ ? Y u + F v u ^ 0 - 


-i av„ 

1 m — 

where ~ -gg- is the derivative in the direction of u. Equa- 
tion (6) becomes 


bn 


h - (v 2 /2) = 2v go cos cp - cos cp 


u 




(7) 
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Because the stagnation enthalpy is assumed to be uniform at the 
entrance to the impeller, equation (2) gives 

Sh = _ 

dn dn 

so that equation (7) becomes 


c3v 

dn 


v v_^ 

^ u m 

2 _ a) cos cp + — 



vr 


cos cp 


1 

v Sn 


If an infinite number of blades is assumed with loading and volume 
dispersed in the fluid, the variation of mass flow is 


- dW = pv ( 2nr — 

m V cos a 


dn 


( 8 ) 


Then 


S(v 2 / 2 ) 

- 

v 




P fenr - Bd ) 

[_ \ cos q/J 


cos q /2 (jl) sin o' 

sin q cos cp 


d 

E m \T r y 

cos q 


dW 


(“l'Ti) 


(9) 


The quantity d(u 1 *V^)/dW is a known function of the particular 
streamline in application of this equation to impeller design. 


The angle \J ; between the axial direction and the curve formed 
by the intersection of the blade with the circular cylinder about 
the axis is related to a and cp by 


tan \|/ 


sin a 

cos cp cos a 


tan a 
cos cp 


If the blades consist of radial elements, 


( 10 ) 


tan \|/ = — tan 

Y» C 

c 
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a 

B 

a 

F 

f 

E 

M 

M* 

n 

P 

R 

r 

S 

3 

T 

u 

V 

V 

W 


APPENDIX B 
SYMBOIS 

The following symbols are used, in the text and the appendix; 
sonic velocity 
number of blades 
thickness of blade 

force exerted on fluid by blades and viscous forces 
slip factor for discharge tangential velocity 
enthalpy of gas 

Mach number of absolute velocity 
Mach number of relative velocity 

coordinate on curve in meridional plane and normal to the 
relative velocity 

pressure 

gas constant 

radius of curvature of meridional projection of streamline 
radial distance from axis of rotation 
entropy of gas 

coordinate along a streamline 
absolute temperature 
linear velocity of rotating impeller 
absolute velocity of gas 
relative velocity of gas 

mas 3 flow between two stream surfaces of revolution 
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a angle "between meridional component and resultant relative gas 
velocity 

7 ratio of specific heats of gas 

d ( r I v u,l )2 

^ aw^ — 

p gas density 

cp angle between meridional and axial components of gas velocity 
'l' tan' 1 (v u /v a ) 

(jl) angular velocity of impeller 
Unit vectors : 

i in direction of relative velocity v 

ft in direction of velocity of impeller u 

ft in meridional plane normal to stream surface of revolution 



m parallel to projection of relative velocity curve on 
meridional plane. (Normal to ft and ft) 

Subscripts : 

1 impeller entrance 

2 inducer exit 

3 impeller discharge 

a axial velocity component 

c case value (at blade tips) 

h hub value (at blade roots) 

Z lift 
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m meridional velocity component 
r radial velocity component 
t stagnation state 

u rotational velocity component 
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ure 2. - Mass flow, pressure, and enthalpy 
at impeller exit. 



Mass flow, Ib/sec Absolute velocity, ft/sec 
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2 3 4 5 

Entrance radius, in. 

Figure 3. - Mass flow and absolute velocity 
at impeller entrance. 
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Abstract 


An axial -discharge centrifugal compressor, 
especially adapted for Jet engines because of the 
large mass flow per unit frontal area, is described. 
General equations of relative fluid motion are 
developed -to show assumptions involved and empirical 
character of simplifications. 

The best impeller was selected on basis of 
maximum air-flow capacity, which was 19.6 pounds 
per second for 14 -inch diameter with tip speed of 
1480 feet per second and pressure ratio of 3.5. 




